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The pathogenesis of severe acute respiratory syndrome-associated coronavirus (SARS-CoV) at the cellular 
level is unclear. No human cell line was previously known to be susceptible to both SARS-CoV and other human 
coronaviruses. Huh7 cells were found to be susceptible to both SARS-CoV, associated with SARS, and human 
coronavirus 229E (HCoV-229E), usually associated with the common cold. Highly lytic and productive rates of 
infections within 48 h of inoculation were reproducible with both viruses. The early transcriptional profiles of 
host cell response to both types of infection at 2 and 4 h postinoculation were determined by using the 
Affymetrix HG-U133A microarray (about 22,000 genes). Much more perturbation of cellular gene transcription 
was observed after infection by SARS-CoV than after infection by HCoV-229E. Besides the upregulation of 
genes associated with apoptosis, which was exactly opposite to the previously reported effect of SARS-CoV in 
a colonic carcinoma cell line, genes related to inflammation, stress response, and procoagulation were also 
upregulated. These findings were confirmed by semiquantitative reverse transcription-PCR, reverse transcrip¬ 
tion-quantitative PCR for mRNA of genes, and immunoassays for some encoded proteins. These transcrip- 
tomal changes are compatible with the histological changes of pulmonary vasculitis and microvascular 
thrombosis in addition to the diffuse alveolar damage involving the pneumocytes. 


Severe acute respiratory syndrome (SARS)-associated coro¬ 
navirus (SARS-CoV) is the etiological agent of SARS (26, 39). 
The disease is associated with significant mortality and mor¬ 
bidity (38). Such aggressive clinical behavior is very different 
from that of other known human coronaviruses, such as the 
group 1 coronaviruses 229E and NL63 or the group 2 corona¬ 
virus OC43. Although these viruses are generally associated 
with mild upper respiratory tract infection such as the common 
cold (47-49), these human coronaviruses can also cause pneu¬ 
monia when the very young, elderly, or immunosuppressed 
hosts are affected (14, 15, 40). Animal models of SARS were 
established by infecting cynomolgus monkeys, ferrets, domestic 
cats, or mice (16, 27, 31, 33, 55), but there are questions about 
their reproducibility of the pathology mimicking human dis¬ 
ease. Although the cellular receptors for attachment of SARS- 
CoV were found to be ACE2 (30) and, recently, L-SIGN (23), 
pathogenesis at the cellular level is largely unknown. 

SARS-CoV is an enveloped, positive-sense, single-stranded 
RNA virus which can grow in embryonal monkey cell lines, 
including Vero E6 and fetal rhesus monkey kidney (FRhk-4) 
cells (26, 39). It can be subcultured onto other Vero cells and 
colonic carcinoma cell lines such as Caco-2 or LoVo (3, 10). 
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Unlike other human coronaviruses, SARS-CoV proliferates 
rapidly and causes obvious cytopathic effects in Vero E6 cells 
within 48 h of inoculation (35). There are no other human cell 
lines known to be susceptible to infection by both SARS-CoV 
and other human coronaviruses. Recently there have been 
reports of a human hepatoma cell line, Huh7, which can be 
infected by pseudotyped lentiviral particles carrying the spike 
protein of SARS-CoV and wild-type replicative SARS-CoV (6, 
20, 45). We report in this study the susceptibility of the cell line 
Huh7 to infection by both SARS-CoV and human coronavirus 
229E (HCoV-229E). A comparative gene transcriptional pro¬ 
file at an early stage of infection of Huh7 cells by these two 
viruses was performed to elucidate their difference and its 
importance in the pathogenesis of disease. 

MATERIALS AND METHODS 

Cell lines and virus. The human hepatoma cell line Huh7 (courtesy of David 
Ho, Aaron Diamond AIDS Research Center) was used throughout this study. 
The cells were incubated at 37°C in minimal essential medium supplemented 
with 10% fetal calf serum and 100 IU/ml of penicillin and 100 |Jig/ml of strep¬ 
tomycin. Our prototype virus (SARS-CoV HKU-39849) was isolated from the 
lung tissue biopsy of the brother-in-law of the index patient who traveled from 
Guangzhou and started a superspreading event in the Hong Kong Special Ad¬ 
ministrative Region leading to the pandemic (37). The HCoV-229E strain 
(ATCC VR-740) was used in this study. The SARS-CoV and HCoV-229E strains 
used in our experiments had undergone three passages in the FRhk-4 cell and 
MRC-5 cell lines, respectively, and were stored at —70°C. Viral titers were 
determined as median tissue culture infective dose (TCID 50 ) per ml in confluent 
Huh7 cells in 96-well microtiter plates, which standardized the viral inoculum 
and measured the relative susceptibility of the Huh7 cell line to these two viruses. 
The relative susceptibilities of the Vero 1008, Vero 76, Vero, and Huh7 cell lines 
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to SARS-CoV and HCoV-229E were also tested by TCID 50 . One hundred 
TCID 50 was confirmed by plaque assays to be equivalent to 85 PFU. All work 
with infectious virus was performed inside a type II biosafety cabinet in a 
biosafety containment level III facility, and the personnel wore powered air- 
purifying respirators (HEPA Airmate; 3 M, Saint Paul, Minn.). 

Monitoring of virus-induced cytopathic effect, antigen detection, and semi- 
quantitative and quantitative RT-PCR. Huh7 cells and culture supernatants 
infected with either SARS-CoV or HCoV-229E at a multiplicity of infection of 
100 TCID 50 per cell were collected at 2, 4, 12, and 24 h postinfection. A washing 
step was performed 1 h postinoculation. The percentages of cells developing 
cytopathic effects were counted by inverted light microscopy at 24 and 48 h. The 
rate of viral replication was measured by reverse transcription-quantitative PCR 
(RT-qPCR) on the culture filtrate. The amount of coronavirus antigen expres¬ 
sion in infected cells was measured by indirect immunofluorescence with conva¬ 
lescent-phase sera of patients with SARS-CoV or HCoV-229E infection, as 
reported previously (2, 8). Briefly, harvested cells were prepared and fixed in 
ice-cold acetone for 10 min. Convalescent-phase serum at a dilution of 1 in 100 
was used to react with the infected cells harvested at various time points. After 
30 min incubation, the cells were washed twice in phosphate-buffered saline for 
5 min each, and then goat anti-human fluorescein isothiocyanate conjugate 
(INOVA Diagnostics, Inc., San Diego, CA) was added and the cells were further 
incubated for 30 min at 37°C. The cells were washed again as described above, 
and the percentage of positive cells was manually estimated under UV micros¬ 
copy. 

Reverse transcription-PCR (RT-PCR) for SARS-CoV and HCoV-229E was 
done directly on culture filtrate according to our previous protocol (2). Briefly, 
total RNA extracted from culture filtrate with the QLAamp virus RNA minikit 
(Qiagen) as instructed by the manufacturer was reverse transcribed with random 
hexamers. cDNA was amplified with SARS-CoV primers (forward, 5'-TACAC 
ACCTCAGCGTTG-3'; reverse, 5'-CACGAACGTGACGAAT-3') and HCoV- 
229E primers (forward, 5'-GGTACTCCTAAGCCTTCTTCG-3'; reverse, 5'-G 
ACTATCAAACAGCATAGCAGC-3'). Using real-time RT-qPCR assays, 
cDNA was amplified in SYBR Green I fluorescence reactions (Roche, Mann¬ 
heim, Germany). For the RT-qPCR of SARS-CoV, a 20-|xl reaction mixtures 
containing 2 p,l cDNA, 3.5 mmol/liter magnesium chloride, and 0.25 |xmol/liter 
of the same forward and reverse primers as in the reaction mixtures were thermal 
cycled with a Light Cycler (Roche) (95°C for 10 min, followed by 50 cycles of 
95°C for 10 s, 57°C for 5 s, and 72°C for 9 s with ramp rates of 20°C/s). For the 
RT-qPCR of HCoV-229E, the conditions were similar to those described above 
except that 45 cycles of 95°C for 10 s, 65°C for 3 s, and 72°C for 12 s were used. 
Plasmids with the target sequences were used to generate the standard curve. At 
the end of the assay, the PCR products of SARS-CoV and HCoV-229E (182 bp 
and 295 bp, respectively) were subjected to a melting curve analysis (65 to 95°C, 
0.1°C/s) to determine the specificity of the assay. All assays were performed in 
replicates. 

Microarray analysis. Human genome-wide gene expression was examined 
with the GeneChip system HG-U133A microarray (Affymetrix Inc., Santa Clara, 
CA), which is composed of more than 22,000 oligonucleotide probe sets inter¬ 
rogating approximately 18,400 unique transcripts, including 14,500 well-charac¬ 
terized human genes. Quality control, GeneChip hybridization, and data acqui¬ 
sition and analysis were performed at the Genome Research Centre, The 
University of Hong Kong, according to the standard protocols available from 
Affymetrix. In brief, total RNAs of the infected or uninfected cell lines at 
different time points were extracted using the RNeasy minikit (Qiagen, Valencia, 
CA). Double-stranded cDNA was synthesized from 10 |jug of total RNA with the 
GeneChipT7-01igo (dT) Promoter Primer Kit (Affymetrix, Inc.) and the Super- 
Script Choice System (Invitrogen). Biotin-labeled cRNA was then synthesized by 
in vitro transcription using the Bio Array High Yield RNA Transcript Labeling 
Kit (Affymetrix, Inc). After fragmentation, 15 |xg of labeled cRNA was hybrid¬ 
ized to the oligonucleotide microarray. The chips were washed and stained using 
the GeneChip Fluidics Station 400 (Affymetrix) and then scanned with the 
GeneChip Scanner 3000 (Affymetrix). Data analysis was performed using the 
Microarray Suite Expression Analysis software (version 5.1; Affymetrix). For 
comparison across different arrays, the data for each array were normalized by a 
global scaling strategy, using a scaling target intensity of 500. By using the 
Affymetrix-defined comparison mathematical algorithms, a fold change in ex¬ 
pression between each of the infected samples in comparison to the uninfected 
mock control was calculated, log 2 transformed, and further classified as not 
changed, increased (signal log ratio change P value of <0.005), decreased (signal 
log ratio change P value of >0.995), or marginally increased or decreased. To 
classify a gene as significantly upregulated or downregulated after infection at a 
specific time point, two additional criteria were used: (i) the fold change must be 
greater than or equal to 2 (signal log ratio of >1 if upregulated or <1 if 


downregulated) to be classified as increased or decreased, and (ii) genes that 
were classified as upregulated must be flagged as present in the infected samples, 
while genes that were classified as downregulated must be flagged as present in 
the uninfected control sample. All gene chip procedures were performed in 
replicates. 

Gene expression analysis by semiquantitative PCR, RT-qPCR, and immuno¬ 
assays. Genes with significant transcriptional changes known to be associated 
with biological significance were selected for further analysis by semiquantitative 
PCR, RT-qPCR, and immunoassays. RT-qPCR was performed according to our 
previous protocol (56). The extracted RNA was pretreated with DNase. Primers 
which specifically amplified nine genes related to apoptosis, inflammation, and 
coagulation were designed (Table 1). First-strand cDNA was synthesized from 
the total RNA by reverse transcription with random hexamers. Semiqualitative 
comparison was performed using simple gel electrophoresis and ethidium bro¬ 
mide staining (10). Quantitative PCR was performed using the SYBR Green I 
fluorescence reactions in a Light Cycler as described above. Detailed PCR 
conditions are available upon request. Serial dilutions of a reference cDNA 
derived from the SARS-CoV-infected sample were used to generate the standard 
curve. Melting curve analysis was performed for each primer pair at the end of 
the reaction to confirm the specificity of the assay. The housekeeping porpho¬ 
bilinogen deaminase gene was used for standardization of the initial RNA con¬ 
tent of a sample. Experiments were performed in duplicate, and the result for an 
individual sample was expressed as the mean expression level of a specific 
gene/porphobilinogen deaminase gene relative to the reference cDNA. The 
relative expression between each infected sample and the uninfected control was 
then calculated and expressed as fold change. 

Three sets of immunoassays (human interleukin-8 [IL-8] from BD Biosciences 
and PAH [serpine 1] and TFPI2 from Diagnostica Stago) were performed 
according to our previous protocols (5, 28) and the manufacturers’ instructions. 
Briefly, for IL-8 assay, 100 jxl culture supernatant or standard was added to IL-8 
monoclonal antibody-coated wells and incubated for 2 h at room temperature. 
After washing five times, 100 |xl a biotinylated anti-human IL-8 monoclonal 
antibody was added to each well. After five washings, 100 p.1 of tetramethyl 
benzidine substrate was added, followed by incubation for 30 min at room 
temperature. The reaction was stopped by adding 50 |xl stop solution. The 
absorbance was read at 450 nm within 30 min. The experimental conditions and 
procedure were similar for the immunoassays for PAI1 and TFP1. All experi¬ 
ments were done in replicates to ensure accuracy and reproducibility. 

Statistical analysis. The fold change of the target gene expressions and the 
difference in the concentrations of expressed proteins between SARS-CoV and 
HCoV-229E at different postinoculation time points were compared by Student’s 
t test. A P value of <0.05 was considered significant. A statistical package (SPSS 
10.0; SPSS Hong Kong, Hong Kong) was used for all analyses. 

RESULTS 

Susceptibility of the Huh7 cell line to SARS-CoV and HCoV- 
229E. With a multiplicity of infection of 100, a cytopathic effect 
was visible in Huh7 cells at 24 and progressed to about 50% 
cell death at 48 h with both viruses. SARS-CoV but not HCoV- 
229E grew well in the Vero cell line. The viruses produced 
comparable TCID 50 of around 10 7 per ml in the culture su¬ 
pernatant of Huh7 cells at 48 h (Table 2). In terms of viral 
load, a 1-log-unit increase of viral genome copy was noted at 
12 h postinoculation with both viruses, which was followed by 
a peak at 24 h (Table 3). Antigen expression could be observed 
by indirect immunofluorescence in over 50% of the cells at 
24 h postinoculation for both viruses. 

Effects on gene expression of host cells determined by mi¬ 
croarray. Two hundred twenty-four genes were significantly 
altered within 4 h postinoculation in the transcriptomal expres¬ 
sion analysis. Only 21 genes were perturbed by HCoV-229E 
per se, whereas 164 genes were altered by SARS-CoV infec¬ 
tion only, and the remaining 39 genes were altered by both 
coronaviruses. Out of the 164 genes with altered expression in 
SARS-CoV, 38 were upregulated and only one was downregu¬ 
lated at both 2 and 4 h postinoculation. Forty-three were up¬ 
regulated and 16 were downregulated at 2 h postinoculation. 
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TABLE 1. Primers used for checking the expression of nine selected genes and a house keeping gene by semiquantitative RT- PCR 

and RT-qPCR 


Gene product 

Primer 

Sequence (5' to 3') 

Length 

(bp) 

PCR product 
size (bp) 

Tissue factor pathway inhibitor 2 

TFPI2-F 

CACCTATACTGGCTGTGGAGG 

21 

236 

TFPI2-R 

CCTCATGCTGTCATATTATTCTTC 

24 


Plasminogen activator inhibitor 1 (serpine 1) 

PAI1-F 

CCGAGGAGATCATCATGGAC 

20 

259 


PAI1-R 

GCCAAGGTCTTGGAGACAGA 

20 


Thrombospondin 1 

THBS1-F 

ATCATGGCTGACTCAGGACC 

20 

229 


THBS1-R 

GGAAGCCAAGGAGAAGTGAT 

20 


Interleukin-8 

IL8-F 

TGGAAGAGAGCTCTGTCTGGA 

21 

215 


IL8-R 

CCAGGAATCTTGTATTGCATCT 

22 


Nuclear factor of kappa light polypeptide gene 

NFKB2-F 

CCTGGCAGGTCTACTGGAGG 

20 

281 

enhancer in B-cells 2 

NFKB2-R 

AAATAGGTGGGGACGCTGT 

19 


Jun B proto-oncogene 

JUNB-F 

CCAGCTCAAACAGAAGGTCATG 

22 

209 

JUNB-R 

GGAGTCCAGTGTGGTTTGCG 

20 


Pleckstrin homology-like domain, family A, 

PHLDA1-F 

CTCTCATCCTCACTCGCACC 

20 

214 

member 1 

PHLDA1-R 

TTGATCCAAGTGAAGACAATAGAA 

24 


Caspase recruitment domain family, member 10 

CARD10-F 

CACGTGGAGGTGACTGAGAA 

20 

274 

CARD10-R 

CTCAGGCCTCACTGCTGCT 

19 


BCL2-associated X protein 

BAX-F 

GATGATTGCCGCCGTGGACA 

20 

175 

BAX-R 

GATGGTTCTGATCAGTTCCGG 

21 


Porphobilinogen deaminase 

PBGD-F 

AAGATGAGAGTGATTCGCGTG 

21 

226 

PBGD-R 

GGTCCACTTCATTCTTCTCCAG 

22 



Forty-nine were upregulated and 17 were downregulated at 4 h 
postinoculation. In contrast, for HCoV-229E infection, only 
one gene was upregulated and no gene was downregulated at 
both 2 and 4 h postinoculation. No gene was upregulated and 
two were downregulated at 2 h postinoculation. Fourteen 
genes were upregulated and four were downregulated at 4 h 
postinoculation. When multiple transcripts of the same gene 
were eliminated and analyzed, genes related to apoptosis (n = 
23), inflammatory or immune response (n = 34), and coagu¬ 
lation (n = 5) were identified, in addition to the expected 
genes related to the stress response and metabolism and other 
unknown genes (Table 4). Of the 23 apoptotic genes affected, 
13 of them are proapoptotic and 11 of them are upregulated in 
SARS-CoV infection, compared with only 3 in HCoV-229E 
infection (Table 5). For inflammation and immune response, 
32 genes are upregulated in SARS-CoV, compared with only 3 
in HCoV-229E. These include the genes for NFKB1A; 


TABLE 2. Viral load of SARS-CoV and HCoV-229E in different 
cell lines 48 h postinoculation 


Cell line 


Log 10 TCID 50 /ml 

SARS Co-V 

HCoV-229E 

Vero 1008 

6.5 

ND" 

Vero 76 

8.33 

ND 

Vero 

6.5 

ND 

Huh7 

7.0 

6.5 


NFKB2; IL-8; transforming growth factor (32; chemokines 
CXCL1, -2, -3, -5, -6, and -10; ICAM1; and tumor necrosis 
factor alpha-induced proteins. Quite unexpectedly, genes of 
the procoagulation pathway were also affected by SARS-CoV 
infection, with upregulation of PLSCR1 (phospholipid scram- 
blase 1), EGR1 (early growth response 1 gene), PAI1/SER- 
PINE1 (plasminogen activator inhibitor 1), and THBS1 
(thrombospondin 1). In terms of stress response, seven genes 
were upregulated in SARS-CoV infection, compared with only 
one in HCoV-229E infection. Overall there were far more 
changes in expression of genes related to the cell cycle, tran¬ 
scription, and metabolism and those with miscellaneous and 
unknown functions in SARS-CoV infection. When the Path¬ 
way Assist software (Ariadne Genomics Inc.) was used for 
linking the altered genes in cellular pathways for SARS-CoV, 
there were clear clusterings of altered genes related to apo¬ 
ptosis, inflammation, and coagulation (Fig. 1). 


TABLE 3. Number of viral genome copies of SARS-CoV and 
HCoV-229E at different postinoculation time points in the Huh7 
cell line 

Copies/ml 


h postinoculation - 

SARS-CoV HCoV-229E 


2 2.2 x 10 7 3.8 X 10 7 

4 1.5 x 10 7 3.3 X 10 7 

12 3.6 x 10 s 7.8 X 10 s 

24 1.2 x 10 9 3.4 X 10 10 


ND, not detectable. 
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TABLE 4. Numbers of genes with significant change of expression 
in different categories and subcategories 


Category and 
subcategory 



No. of genes 


Total 

SARS- 

CoV 

HCoV- 

229E 

Both SARS-CoV 
and HCoV-229E 

Apoptosis 

23 




Antiapoptosis 

3 

3 

0 

0 

Proapoptosis 

12 

7 

2 

3 

DNA repair 

5 

2 

2 

1 

Other 

3 

2 

1 

0 

Inflammation 

34 




Antiinflammatory 

2 

2 

0 

0 

Proinflammatory 

24 

21 

0 

3 

Immunomodulatory 

2 

1 

1 

0 

Other 

6 

5 

1 

0 

Coagulation 

5 




Anticoagulants 

1 

0 

0 

1 

Procoagulants 

4 

4 

0 

0 

Metabolism 

35 




Amino acid 

4 

2 

1 

1 

Carbohydrate 

7 

6 

0 

1 

Lipid 

5 

5 

0 

0 

Transporter 

10 

6 

1 

3 

Other 

9 

7 

2 

0 

Signal Transduction 

18 

15 

1 

2 

Transcription 

37 

25 

3 

9 

Cell Cycle 

17 

12 

3 

2 

Cell growth and 

7 

6 

0 

1 

differentiation 





Stress response 

7 

6 

0 

1 

Miscellaneous 

16 

11 

1 

4 

Unknown 

25 

16 

2 

7 

Total 

224 

164 

21 

39 


Confirmation of cellular gene and protein expression by 
semiquantitative PCR, RT-qPCR, and immunoassays. Since 
previous studies have documented the alterations of the apo- 
ptotic and inflammatory pathways in relation to SARS-CoV 
infection, only key genes related to inflammation and apopto¬ 
sis were chosen to confirm the results of microarray analysis by 
semiquantitative RT-PCR and RT-qPCR. For the coagulation 
pathway, both procoagulation and anticoagulation genes were 
chosen. A similar trend of upregulation of gene expression was 
found by RT-qPCR, with SARS-CoV showing a 1.4- to 10.8- 
fold increase compared with HCoV-229E for coagulation 
(TFPI2, PAI1, and THBS1), inflammation (IL-8 and NFKB2), 
transcription (JUNB), and apoptotic (PHLDA1, CARD 10, 
and BAX) pathways (Fig. 2 and 3). IL-8 was increased in both 
SARS-CoV and HCoV-229E infection, but the increase was 
much higher in the case of SARS-CoV. The enzyme immuno¬ 
assay showed that SARS-CoV induced higher concentrations 
of PAI1 and IL-8 than HCoV-229E at 2, 4, 12, and 24 h 
postinoculation. Both SARS-CoV and HCoV-229E induced 
similar TFPI2 expression at 4,12, and 24 h postinoculation, but 


at 2 h postinoculation SARS-CoV induced a lower concentra¬ 
tion of this predominately anticoagulation protein (Fig. 4). 

DISCUSSION 

Of the four coronaviruses known to infect human, HCoV- 
229E, HCoV-OC43, and NL63 are generally associated with 
mild upper respiratory tract infection such as the common cold 
in immunocompetent hosts (47-49). In contrast, SARS-CoV 
causes respiratory failure in over 60% of affected persons, with 
a mortality rate of 15% (37, 38). Besides pneumonia, SARS is 
also clinically manifested as watery diarrhea without enteroco¬ 
litis (48.6%) (8), hepatitis without liver failure (49.4%) (22), 
lymphopenia (75%) (37), impaired coagulation (63%) (51), 
and occasionally pulmonary vasculitis and thrombosis in the 
lungs of those who died (36, 13). Much is already known about 
many aspects of SARS, including the virology (18, 26, 39, 53), 
genomics (32), diagnostics (5, 28, 43, 52, 54), clinical features 
and progression in relation to viral load (22, 37), treatment (7, 
9, 46), infection control (44), and immunization (55). However 
little is known about the pathogenesis at the cellular level, 
despite the identification of ACE2 and L-SIGN as the recep¬ 
tors for binding of SARS-CoV to host cells (30, 23). 

One important limitation to studies at the cellular level is the 
absence of a relevant cell line which can be lytically infected by 
SARS-CoV. Initially, only embryonal monkey kidney cell lines 
such as Vero or FRhk-4 and their derivatives could be readily 
infected. Subsequent investigation showed that human colon 
adenocarcinoma cell lines Caco-2, CL-14, and LoVo could also 
be infected (3, 10). However, inflammation or cell death is not 
manifested on endoscopic examination of the intestines of 
SARS patients (29). Examination of intestinal tissue biopsy by 
electron microscopy revealed abundant intracellular viral par¬ 
ticles, and there was negligible inflammatory cells or cellular 
apoptosis on light microscopy (29). In fact, a high TCID 50 in 
the culture supernatant of persistently infected LoVo cells can 
be achieved without any cytopathic effects (3). Recently a gene 
expression analysis showed that SARS-CoV upregulated anti- 
apoptotic genes and several CXC chemokines while downregu¬ 
lating proapoptotic genes, IL-18, and macrophage migration 
inhibitory factor (10) at 24 h postinoculation of Caco-2 cell 
lines. These patterns of gene expression appeared to explain 
the severe watery diarrhea without clinical or endoscopic signs 
of necrosis or inflammation in SARS patients (29). However, 
no data on early gene expression before the onset of cytolysis 
were reported. Similarly, no comparative study of the host cell 
transcriptional profile and those for other, less virulent coro¬ 
naviruses was ever reported. 

No pneumocyte cell line has yet been found to support lytic 
or nonlytic infection by SARS-CoV. However, the human hep¬ 
atoma cell line Huh7 was found to be susceptible to SARS- 
CoV. This is not completely unexpected, because the mouse 
hepatitis virus, a group 2 coronavirus, is known to infect Huh7 
cells (25). In this study HCoV-229E was found to produce 
good lytic infection within 48 h postinoculation. A high multi¬ 
plicity of infection of 100 TCID 50 s per cell was used to ensure 
the reproducibility of gene expression, as previously reported 
(41). Since the expression of a high number of genes was 
expected to change significantly when virus-induced cytopa- 
thology was impending for a rapidly lytic viral infection, we 
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TABLE 5. Signal log ratios (log 2 scale) of the genes that are significantly altered after SARS-CoV or HCoV-229E infection 


Category and 
subcategory 


Gene 


Signal log ratio (log 2 scale)" 
after infection with: 

Asymetrix 

identification 

no. 

Designation 

Description 

HCoV-229E 

SARS-CoV 



2 h 

4 h 

2 h 

4 h 

Apoptosis 

Antiapoptosis 

200797 s at 

MCLl 

Myeloid cell leukemia sequence 1 (BCL2 related) 

0.2 

0 

1 

0.3 


210538 s at 

BIRC3 

Baculoviral IAP repeat containing 3 

0.2 

0.6 

2.1 

1.4 


201631_s_at 

IER3 

Immediate-early response 3 

0.3 

-0.2 

L7 

L3 

Apoptosis 

204285 s at 

PMAIP1* 

Phorbol-12-myristate-13-acetate-induced protein 1 

0.5 

-0.2 

1.6 

0.6 


204908 s at 

BCL3 fc 

B-cell CLL/lymphoma 3 

0.1 

0.1 

1.9 

1.1 


207181 s at 

CASP7 

Caspase 7, apoptosis-related cysteine protease 

-0.2 

0.4 

0.1 

1.2 


208478 s at 

BAX 

BCL2-associated X protein 

0.5 

0.6 

0.4 

1.1 


208536 s at 

BCL2L11 

BCL2-like 11 (apoptosis facilitator) 

-1.3 

0.3 

0.1 

1.2 


209305 s at 

GADD45B 6 

Growth arrest and DNA-damage-inducible, beta 

0.1 

-0.2 

1.4 

0.8 


210025 s at 

CARD10 6 

Caspase recruitment domain family, member 10 

1 

-0.4 

1.9 

2.7 


211085 s at 

STK4 

Serine/threonine kinase 4 

0.4 

1.3 

0.2 

0.6 


211695 x at 

MUC1 

Mucin 1. transmembrane 

-0.2 

-1.3 

-0.1 

-0.8 


218368_s_at 

TNFRSF12A 

Tumor necrosis factor receptor superfamily, 
member 12A 

0.5 

0 

i 

0.6 


218880 at 

FOSL2 

FOS-like antigen 2 

1.8 

-0.4 

3.3 

1.5 


217997_at* 

PHLDA1 

Pleckstrin homology-like domain, family A, 
member 1 

0.5 

0.1 

2.9 

L6 

DNA repair 

205022 s at 

CHES1 

Checkpoint suppressor 1 

0 

-0.2 

0.2 

-1 


205072_s_at 

XRCC4 

X-ray repair complementing defective repair in 
Chinese hamster cells 4 

-1.1 

-0.2 

-1.3 

-0.3 


208386_x_at 

DMC1 

DMC1 dosage suppressor of mckl homolog, 
meiosis-specific homologous recombination 

-0.1 

L2 

0.4 

0.2 


209579 s at 

MBD4 

Methyl-CpG binding domain protein 4 

0.1 

0.6 

0.2 

1 


219317_at 

POLI 

Polymerase (DNA directed) iota 

0 

LI 

-0.6 

0.5 

Other 

204178 s at 

RBM14 

RNA binding motif protein 14 

-0.3 

0.5 

-0.3 

1 


209803_s_at 

PHLDA2 

Pleckstrin homology-like domain, family A, 
member 2 

0.1 

0.1 

0.4 

I 


215338_s_at 

NKTR 

Natural killer-tumor recognition sequence 

0.8 

1 

0.3 

0.2 

Inflammation 

Anti-inflammatory 

201502_s_at 

NFKBIA/IKBA 

Nuclear factor of kappa light polypeptide gene 
enhancer in B-cells inhibitor, alpha/I-kappa B 
alpha 

0.2 

-0.1 

L4 

0.5 


221577_x_at 

PLAB 

Prostate differentiation factor 

0.6 

-0.2 

L8 

0.5 

Proinflammatory 

202510 s at 

TNFAIP2 

Tumor necrosis factor alpha-induced protein 2 

-0.3 

-0.3 

1.2 

1.2 


202625_at 

LYN 6 

v-yes-1 Yamaguchi sarcoma viral related 
oncogene homology/oncogene LYN 

0.1 

-0.1 

-0.2 

L2 


202637_s_at 

ICAMl 6 

Intercellular adhesion molecule 1 (CD54), human 
rhinovirus receptor 

0.2 

-0.3 

2.5 

L8 


202644 s at 

TNFAIP3'’ 

Tumor necrosis factor alpha-induced protein 3 

0.2 

-0.1 

2.1 

1 


203313 s at 

TGIF 

TGFB-induced factor (TALE family homeobox) 

0.5 

0.3 

1 

1 


203964_at 

NMI 

N-myc (and STAT) interactor (transcription from 
Pol II promoter, inflammatory response, JAK- 
STAT cascade) 

0 

0.2 

0.2 

I 


204470_at 

CXCL1 

Chemokine (C-X-C motif) ligand 1 (melanoma 
growth-stimulating activity, alpha) 

0.8 

-0.2 

3.1 

0.9 


205114_s_at 

CCL3/MIP1A 

Chemokine (C-C motif) ligand 3/macrophage 
inflammatory protein 1-alpha 

0.5 

0.4 

L4 

1 


205476_at 

CCL20/MIP3A 

Chemokine (C-C motif) ligand 20/macrophage 
inflammatory protein 3-alpha 

0.4 

-0.2 

L2 

0.8 


205844 at 

VNN1 

Vanin 1 

-0.6 

0.3 

0.1 

1.2 


205922 at 

VNN2 

Vanin 2 

0.7 

-0.2 

0.8 

1 


206336_at 

CXCL6 

Chemokine (C-X-C motif) ligand 6 (granulocyte 
chemotactic protein 2) 

0.2 

0.4 

0.5 

L3 


207850_at 

CXCL3/MIP2B 

Chemokine (C-X-C motif) ligand 3/macrophage 
inflammatory protein 2-beta 

0.2 

L6 

5.2 

3.5 


208991_at 

STAT3 

Signal transducer and activator of transcription 3 
(acute-phase response factor) 

0 

-0.5 

1 

0.4 
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TABLE 5 —Continued 


Category and 
subcategory 


Gene 


Signal log ratio (log 2 scale) 0 
after infection with: 

Affymetrix 

identification 

no. 

Designation 

Description 

HCoV-229E 

SARS-CoV 



2 h 

4 h 

2 h 

4 h 


209636_at 

NFKB2 6 

Nuclear factor of kappa light polypeptide gene 
enhancer in B cells 2 

-0.8 

-1.4 

2.5 

2.7 


209774_x_at 

CXCL2/MIP2A 

Chemokine (C-X-C motif) ligand 2/macrophage 
inflammatory protein 2-alpha 

IT 

0 

3 

L6 


209909 s at 

TGFB2 6 

Transforming growth factor beta 2 

0.1 

0.2 

0.9 

2.1 


211000 s_at 

IL-6ST 6 

Interleukin-6 signal transducer (gpl30, oncostatin 

M receptor) 

0.5 

0.9 

0 

L3 


211506 s at 

IL-8 6 

Interleukin-8 

0.6 

0.7 

3.3 

2.1 


212063_at 

CD44 

CD44 antigen (homing function and Indian blood 
group system) 

0.1 

-0.1 

0.5 

IT 


212989_at 

CXCL10 

Chemokine (C-X-C motif) ligand 10 or Mob 
protein 

0.2 

0.6 

0.9 

I 


215101 s at 

CXCL5 

Chemokine (C-X-C motif) ligand 5 

-0.2 

0.4 

1.6 

2 


217738 at 

PBEF 6 

Pre-B-cell colony-enhancing factor 

0 

-0.1 

1.5 

1.3 


AFFX- 

HUMIS 

GF3A/ 

M97935 

_5_at 

ISGF-3 

Interferon-stimulated transcription factor 3 

0.1 

0.4 

0.4 

L3 

Immuno¬ 

modulatory 

201667_at 

GJA1/CX43 

Gap junction protein, alpha 1/gap junction 
protein 43-KD/connexin 43/heart connexin 

0 

0.3 

0.7 

L4 


220954_s_at 

PILR, BETA 

Paired immunoglobulin-like receptor beta 

-0.2 

1 

-0.3 

0.9 

Other 

202531 at 

IRF1 

Interferon regulatory factor 1 

0.1 

-0.6 

1 

0.1 


204622 x at 

NR4A2 

Nuclear receptor subfamily 4, group A, member 2 

0.3 

-0.3 

-1.6 

-0.9 


204897 at 

PTGER4 

Prostaglandin E receptor 4 (subtype EP4) 

0.2 

0.2 

1.2 

1.1 


205220_at 

Chemokine receptor 
HM74 

Putative chemokine receptor 

0.2 

-0.1 

L7 

1 


206785_s_at 

KLRC2 

Killer cell lectin-like receptor subfamily C, 
member 2 

0.3 

2.2 

3.8 

5.2 


22049l_at 

HAMP 

Hepcidin antimicrobial peptide 

0.4 

-0.4 

2.1 

0.2 

Coagulation 

Anticoagulants 

209277_at 

TFPI2* 

Tissue factor pathway inhibitor 2 

1 

0.9 

2.1 

3 

Procoagulants 

202446 s at 

PLSCR1/MMTRA1B 

Phospholipid scramblase 1 

0 

0.2 

0.6 

1.2 


201110 s at 

THBSl 6 

Thrombospondin 1 

0.5 

0.2 

2.3 

1.1 


201693_s_at 

EGR1 6 

Early growth response 1 

0.1 

-0.9 

1 

-1.1 


202627_s_at 

PAI1 

Plasminogen activator inhibitor 1, serpine 1 

-0.4 

-0.6 

1.3 

1.8 

Metabolism 

Amino acid 

201739 at 

SGK 

Serum/glucocorticoid regulated kinase 

0.2 

-0.2 

1 

0.4 


206669 at 

GAD1 

Glutamate decarboxylase 1 (brain, 67 kDa) 

0.1 

2.9 

2.8 

3 


219795_at 

SLC6A14 

Solute carrier family 6 (neurotransmitter 
transporter), member 14 

-0.1 

0.4 

IT 

2 


222018_at 

NACA 

Nascent-polypeptide-associated complex alpha 
polypeptide 

0.3 

1 

-0.4 

0.4 

Carbohydrate 

202464_s_at 

PFKFB3 

6-Phosphofructo-2-kinase/fructose-2,6- 
biphosphatase 3 

0.1 

0.3 

IT 

IT 


214850 at 

SMA5 

H. sapiens beta glucuronidase pseudogene 

1 

1.5 

0.8 

2.1 


215977 x at 

GK 

Glycerol kinase 

0 

0.6 

0.8 

1 


218918 at 

MAN1C1 

Mannosidase, alpha, class 1C, member 1 

-0.3 

0.2 

-0.8 

-1 


218985_at 

SLC2A8 

Solute carrier family 2 (facilitated glucose 
transporter), member 8 

-0.1 

0.4 

0.3 

I 


219508 at 

GCNT3 

Glucosaminyl (/V-acetyl) transferase 3, mucin type 

0.2 

0.1 

1 

1 


219797_at 

MGAT4A 

Mannosyl (alpha-1,3-)-glycoprotein beta-1,47V- 
acetylglucosaminyltransferase, isoenzyme A 

0.4 

0.7 

0.9 

L2 
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Category and 
subcategory 


Gene 


Signal log ratio (log 2 scale) 0 
after infection with: 

Affymetrix 

identification 

no. 

Designation 

Description 

HCoV-229E 

SARS-CoV 



2 h 

4 h 

2 h 

4 h 

Lipid 

202067_s_at 

LDLR* 

Low-density lipoprotein receptor (familial 
hypercholesterolemia) 

0.5 

-0.3 

L5 

-0.3 


204242 s at 

ACOX3 

Acyl coenzyme A oxidase 3, pristanoyl 

0.3 

0 

0.5 

1.3 


205769_at 

SLC27A2 

Solute carrier family 27 (fatty acid transporter), 
member 2 

0 

-0.3 

-0.6 

-1.2 


205822_s_at 

HMGCS1 

3-Hydroxy-3-methylglutaryl coenzyme A synthase 

1 (soluble) 

0 

0.8 

0.7 

I 


212226_s_at 

PPAP2B 6 

Phosphatidic acid phosphatase type 2B 

0 

-0.3 

L2 

-0.3 

Transporters 

202437_s_at 

CYP1B1 6 

Cytochrome P450, family 1, subfamily B, 
polypeptide 1 

0.6 

L5 

3 

3.5 


204685 s at 

ATP2B2 

ATPase, Ca 2+ transporting, plasma membrane 2 

0 

-0.3 

0.1 

-1 


205195_at 

AP1S1* 

Adaptor-related protein complex 1, sigma 1 
subunit 

L7 

L7 

i 

2.1 


205749_at 

CYP1A1 

Cytochrome P450, family 1, subfamily A, 
polypeptide 1 

-0.4 

L7 

L6 

2.5 


207604_s_at 

SLC4A7 6 

Solute carrier family 4, sodium bicarbonate 
cotransporter, member 7 

0.3 

0.6 

0.7 

L4 


208914_at 

GGA2 

Golgi-associated, gamma adaptin ear-containing, 
ARF binding protein 2 

-0.2 

-0.5 

-0.2 

zY 


210357 s at 

SMOX 

Spermine oxidase 

0 

-0.2 

1.2 

1.3 


211572_s_at 

SLC23A1 

Solute carrier family 23 (nucleobase 
transporters), member 1 

0.3 

0.1 

0 

I 


218703 at 

SEC22L2 

SEC22 vesicle trafficking protein-like 2 

1 

1.3 

-0.5 

0.9 


220786_s_at 

SLC38A4 

Solute carrier family 38, member 4 

-0.2 

0 

-0.7 

-1.1 

Miscellaneous 

201556_s_at 

VAMP2 

Vesicle-associated membrane protein 2 
(synaptobrevin 2) 

0 

-1.8 

-0.7 

-0.8 


202238 s at 

NNMT* 

Nicotinamide V-methyltransferase 

0.4 

0.5 

1.9 

2.8 


203798 s at 

VSNL1 

Visinin-like 1 

-0.1 

-0.1 

0 

-1 


206293_at 

SULT2A1 

Sulfotransferase family, cytosolic, 2A, 

dehydroepiandrosterone preferring, member 1 

-0.2 

0 

-0.3 

:zl 


208776_at 

PSMD11 

Proteasome (prosome, macropain) 26S subunit, 
non-ATPase, 11 

0 

-0.3 

-0.2 

zY 


210447 at 

GLUD2 

Glutamate dehydrogenase 2 

-0.2 

-0.1 

-0.9 

-1 


213926 s at 

HRB 

HIV-1 C Rev binding protein 

0.2 

0.4 

0.3 

1 


214837 at 

ALB 

Albumin 

-0.2 

1.2 

-0.7 

0.2 


216049_at 

RHOBTB3 

Rho-related BTB domain-containing 3 

-0.5 

-0.4 

- 1 

-0.6 

Signal transduction 

201010 s at 

TXNIP 

Thioredoxin-interacting protein 

-0.5 

-0.2 

-1.6 

-0.6 


201044 x at 

DUSP1* 

Dual-specificity phosphatase 1 

0.1 

-0.4 

3.2 

2.1 


201295 s at 

WSB1 

SOCS box-containing WD protein SWiP-1 

0.2 

1.6 

-0.2 

1.4 


201824 at 

RNF14 

Ring finger protein 14 

-0.2 

-0.1 

-1.2 

-0.3 


203879_at 

PIK3CD 

Phosphoinositide-3-kinase, catalytic, delta 
polypeptide 

-0.2 

-0.3 

0.1 

L2 


205205_at 

RELB 

v-Rel reticuloendotheliosis viral oncogene 
homolog B, nuclear factor 
of kappa light polypeptide gene enhancer in B 
cells 3 (avian) 

0.3 

-0.5 

L4 

0.8 


205239_at 

AREG 

Amphiregulin (schwannoma-derived growth 
factor) 

0.4 

0.8 

L5 

0.7 


205868_s_at 

PTPN11 

Protein tyrosine phosphatase, non-receptor type 

11 (Noonan syndrome 1) 

0.6 

L6 

0.6 

0.8 


210056 at 

RH06 

GTP binding protein 

0.7 

0.3 

3.8 

2.6 


211962 s at 

ZFP36L1 

Zinc finger protein 36 

0.2 

0 

1.1 

0.7 


212099 at 

ARHB 

Ras homolog gene family, member B 

0.2 

-0.5 

1.3 

0.6 


212777 at 

SOS1 6 

Son of sevenless homolog 1 (Drosophila) 

0 

0.8 

0.6 

1.2 


213848 at 

DUSP7 

Dual-specificity phosphatase 7 

-0.3 

0.1 

-1 

-0.2 


213880 at 

GPR49 

G protein-coupled receptor 49 

0.2 

0.5 

-0.4 

1.7 


218295 s at 

NUP50 

Nucleoporin, 50 kDa 

0.2 

0.8 

0.4 

1 


219677 at 

SSB1 

SPRY domain-containing SOCS box protein SSB-1 

0.4 

0.2 

1.4 

0.6 
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TABLE 5— Continued 


Gene 

Category and - 

subcategory Aflymetrix 

identification Designation 

no. 


Signal log ratio (log 2 scale) 0 
after infection with: 

Description HCoV-229E SARS-CoV 


2 h 4 h 2 h 4 h 



41660 at 

CELSR1 

Cadherin epidermal growth factor LAG seven- 

0.7 

0.6 

0.4 

1.1 




pass G-type receptor 1 






AFFX- 

GPR34 

G protein-coupled receptor 34 

-2.7 

3 

0.9 

3.7 


HUM 

RGE/ 








M10098_ 

5_at 







Transcription 

201169 s at 

BHLHB2 6 

Basic helix-loop-helix domain containing, class B, 2 

0 

0 

1.3 

1.8 


201473 at 

JUNB 

Jun B proto-oncogene 

0.5 

-0.3 

2.4 

2.2 


201510 at 6 

ELF3* 

E74-like factor 3 (Ets domain transcription 

0.5 

0 

1.9 

1.3 




factor, epithelial specific) 






202601 s at 

HTATSF1/TATSF1 

HIV TAT specific factor 1 

0.7 

1.4 

0.9 

1.1 


202672 s at 

ATF3 

Activating transcription factor 3 

0.6 

-0.6 

2.2 

0.1 


202935 s at 

SOX9 

SRY (sex determining region Y)-box 9 

0.4 

-0.3 

1.1 

-0.3 




(campomelic dysplasia, autosomal sex reversal) 






203751 x at 

JUND 6 

Jun D proto-oncogene 

0.5 

-1.2 

1.4 

-0.4 


203973 s at 

CEBPD 

CCAAT/enhancer binding protein (C/EBP), delta 

0.6 

-0.6 

2.3 

1.5 


204069 at 

MEIS1 

Meisl, myeloid ecotropic viral integration site 1 

-1.1 

0.3 

-0.7 

0.1 




homolog (mouse) 






205187 at 

SMAD5 

Mothers against DPP homolog 5 (Drosophila) 

0.9 

1 

1.3 

1.4 




SMAD 






205398 s at 

SMAD3 

SMAD, mothers against DPP homolog 3 

0.4 

0.1 

1 

1.1 




( Drosophila ) 






205443 at 

SNAPC1 

Small nuclear RNA activating complex, 

0.4 

0.6 

1 

0.8 




polypeptide 1 






205471 s at 

DACH1 

Dachshund homolog (Drosophila) 1 

0.4 

-0.6 

-0.1 

-1 


205975 s at 

HOXD1 

Homeo box D1 

0 

-0.3 

1 

0.7 


206036 s at 

REL 

v-Rel reticuloendotheliosis viral oncogene 

0.4 

0.2 

1 

0.1 




homolog (avian) 






206175 x at 

ZNF222 

Zinc finger protein 222 

0.2 

1 

0.1 

0.5 


207147 at 

DLX2 

Distal-less homeo box 2 

0.7 

-1.8 

1.5 

1.7 


207558 s at 

PITX2 

Paired-like homeodomain transcription factor 2 

-0.2 

1.1 

-1 

-0.1 


208129 x at 

RUNX1 

Runt-related transcription factor 1 (acute myeloid 

0.5 

0.7 

1 

0.4 




leukemia 1; amll oncogene) 






209291 at 

ID4 

Inhibitor of DNA binding 4. dominant negative 

0.2 

-0.2 

-0.6 

-1 




helix-loop-helix protein 






209357 at 6 

CITED2 fc 

Cbp/p300-interacting transactivator, with Glu/ 

0.1 

-0.1 

-1.4 

-0.2 




Asp-rich carboxy-terminal domain, 2 






209505 at 

NR2F1 

Nuclear receptor subfamily 2, group F, member 1 

-0.1 

-0.2 

-1 

-0.3 


209757_s_at 

MYCN 

v-Myc myelocytomatosis viral related oncogene, 
neuroblastoma 
derived (avian) 

0.2 

-0.3 

-1.1 

-0.7 


212641 at 

HIVEP2 

HIV-2 1 enhancer-binding protein 2 

0.3 

-0.1 

1.8 

0.2 


213668 s at 

SOX4 6 

SRY (sex-determining region Y) box 4 

0.5 

0.7 

1.8 

1.3 


213844 at 

HOXA5 

Homeo box A5 

-0.1 

0.1 

1.1 

1 


213906 at 

MYBL1 

v-Myb myeloblastosis viral oncogene homolog 

-0.5 

0.4 

-1 

0.2 




(avian)-like 1 






214639 s at 

HOXA1 

Homeo box A1 

-1.1 

-1.2 

1.3 

1.9 


215073 s at 

NR2F2* 

Nuclear receptor subfamily 2, group F, member 2 

0.1 

0 

-1.2 

-0.5 


215720 s at 

NFYA 

Nuclear transcription factor Y, alpha 

0.4 

1.6 

0.8 

1.4 


218559 s at 

MAFB 

v-Maf musculoaponeurotic fibrosarcoma 

-0.4 

0 

-1.4 

-1.3 




oncogene homolog B (avian) 






220444 at 

ZNF557 

Zinc finger protein 557 

2.3 

2.2 

0.9 

2.4 


221773 at 

ELK3 

ETS domain protein 

-0.3 

0.1 

-0.1 

1 


221883 at 

PKNOX1 

PBX/knotted 1 homeobox 1 

0.8 

0.9 

0.2 

1.7 


36711 at 

MAFF 6 

v-Maf musculoaponeurotic fibrosarcoma (avian) 

0.6 

-1 

1.9 

0 




oncogene family 






40837 at 

TLE2 

Transducin-like enhancer of split 2 [E(spl) 

-0.6 

-1 

-0.5 

-0.5 




homolog, Drosophila] 






44783_s_at 

HEY1 

Hairy/enhancer-of-split related with YRPW motif 1 

0 

0.3 

0.6 

L3 
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Category and 
subcategory 


Gene 


Signal log ratio (log 2 scale) 0 
after infection with: 

Affymetrix 

identification 

no. 

Designation 

Description 

HCoV-229E 

SARS-CoV 



2 h 

4 h 

2 h 

4 h 

Cell cycle 

201055 s at 

HNRPA0 

Heterogeneous nuclear ribonucleoprotein A0 

0.1 

0.6 

0.5 

1 


202150_s_at 

NEDD9 

Neural precursor cell expressed, developmental^ 
downregulated 9 

0 

-0.8 

IT 

0.8 


202768_at 

FOSB 

FBJ murine osteosarcoma viral oncogene 
homolog B 

-0.3 

-0.4 

3.8 

-0.1 


202769 at 

CCNG2 6 

Cyclin G2 

0.1 

0 

-0.2 

1.5 


203120 at 

TP53BP2 

Tumor protein p53 binding protein, 2 

0.3 

0.2 

1.1 

0.4 


204137_at 

TM7SF1 

Transmembrane 7 superfamily member 1 
(upregulated in kidney) 

L7 

L6 

L3 

2.7 


204159_at 

CDKN2C 

Cyclin-dependent kinase inhibitor 2C (pl8, 
inhibits CDK4) 

-0.2 

0.1 

-1.3 

-0.7 


205027 s at 

MAP3K8 

Mitogen-activated protein kinase kinase kinase 8 

-0.1 

0.7 

1.9 

0.7 


205046 at 

CENPE 

Centromere protein E 

0.1 

1.1 

-0.1 

0.8 


205425 at 

HIP1 

Huntingtin-interacting protein 1 

-0.1 

-1 

-0.7 

-0.9 


207530_s_at 

CDKN2B 

Cyclin-dependent kinase inhibitor 2B (pl5, 
inhibits CDK4) 

0.4 

2.3 

2.5 

3.9 


209680 s at 

KIFC1 

Kinesin family member Cl 

0 

0.4 

-0.6 

1 


210136 at 

MBP 

Myelin basic protein 

0.4 

0.5 

1.1 

0.9 


212000 at 

SFRS14 

Arginine/serine-rich 14 splicing factor 

0.8 

1.2 

0.4 

0.8 


212769_at 

TLE3 

Transducin-like enhancer of split 3 [E(spl) 
homolog, Drosophila] 

0 

0.2 

L2 

0.7 


221841 s at 

KLF4 

Kruppel-like factor 4 

0.4 

-0.8 

1.3 

0 


222303_at 

ETS2 

v-Ets avian erythroblastosis virus E26 oncogene 
homolog 2 

-0.2 

0.5 

12 

0.8 

Cell growth and 
differentiation 

201329_s_at 

ETS2 fc 

v-Ets erythroblastosis virus E26 oncogene 
homolog 2 (avian) 

0.3 

0.1 

2 

L3 


201341_at 6 

ENC1/PIG10/NRPB 

Ectodermal-neural cortex (with BTB-like domain) 
p53-induced gene 10; nuclear restricted 
protein/brain 

0 

0.2 

0.8 

L2 


207401 at 

PROX1 

Prospero-related homeobox 1 

0.3 

0.5 

1.2 

1.1 


210001 s at 

SOCS1 

Suppressor of cytokine signaling 1 

0.2 

-0.5 

1.1 

0.4 


210587 at 

INHBE 

Inhibin, beta E 

0.2 

-0.1 

0.9 

1.1 


213004 at 

ANGPTL2 

Angiopoietin-like 2 

-0.1 

0.1 

0.3 

1 


216684_s_at 

SS18 fc 

Synovial sarcoma translocation, chromosome 18 

0.9 

L5 

L3 

L7 

Stress response 

200664 s at 

DNAJB1 

Heat shock 40-kDa protein 1 

0.1 

0.4 

1 

0.7 


202376_at 

AACT/SERPINA3 

Alpha-1 antichymotrpysin/serine (or cysteine) 
proteinase inhibitor, clade A (alpha-1 
antiproteinase, antitrypsin), member 3 

-0.4 

-0.8 

IT 

L3 


20258l_at 

HSPA1B/HSP70-2 

Heat shock 70-kDa protein lB/heat shock 
protein, 70 kDa 2 

-0.2 

0.8 

-0.2 

I 


208607 s at 

SAA2 

Serum amyloid A2 

2.7 

2 

3.5 

3.7 


209799_at 

PRKAA1 

Protein kinase, AMP activated, alpha 1 catalytic 
subunit 

0 

0.7 

0.2 

IT 


215078 at 

SOD2 6 

Superoxide dismutase 2, mitochondrial 

0.4 

0.4 

3 

1.7 


221477_s_at 

SOD7 

Superoxide dismutase 2, mitochondrial 

0 

0.1 

IT 

L3 

Miscellaneous 

202241 at 

TR1B1 

Tribbles homolog 1 ( Drosophila ) 

-0.1 

-0.6 

1.1 

0.3 


202458 at 

SPUVE 

Protease, serine, 23 

0.1 

0.2 

0.4 

1 


202760 s at 

AKAP2 

A kinase (PRKA) anchor protein 2 

0.6 

0 

1.1 

0.1 


203234 at 

UPP1 

Uridine phosphorylase 1 

-0.6 

0.3 

0.3 

1.3 


205890 s at 

UBD 

Ubiquitin D 

0 

-0.2 

1.3 

2 


206374 at 

DUSP8 

Dual-specificity phosphatase 8 

0.6 

-0.6 

1.7 

-0.7 


210159 s at 

TRIM31* 

Tripartite motif-containing 31 

2.2 

1.1 

2.3 

3 


212027 at 

RBM25 

RNA binding motif protein 25 

0.2 

1.1 

0.3 

0.9 


212244_at 

GRINL1A 

Glutamate receptor, ionotropic, IV-methyl D- 
aspartate-like 1A 

-0.2 

-0.1 

-1 

-0.2 


212905_at 

CSTF2T 

Cleavage stimulation factor, 3' pre-RNA, subunit 

2, 64 kDa, tau variant 

-0.5 

0 

-1.1 

0.1 


Continued on facing page 
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TABLE 5 —Continued 


Category and 
subcategory 


Gene 

Description 

Signal log ratio (log 2 scale) 0 
after infection with: 

Affymetrix 

identification 

no. 

Designation 

HCoV-229E 

SARS-CoV 

2 h 

4 h 

2 h 

4 h 


213006 at 

KIAA0146 

KIAA0146 protein 

0.2 

-0.6 

1.8 

1.1 


217475 s at 

LIMK2 

LIM domain kinase 2 

0.2 

0.3 

1 

0.5 


218737 at 

SBNOl 

Sno, strawberry notch homolog 1 

0.9 

1 

0.3 

1 


218810 at 

FLJ23231 

Hypothetical protein FLJ23231 

-1 

-3.5 

2.7 

1.3 


220992 s at 

Clorf25 

Chromosome 1 open reading frame 25 

0.1 

1.2 

-0.1 

1 


221258_s_at 

DKFZP434G2226 

Kinesin family member 18A 

0 

0.6 

-0.9 

I 

Unknown 

204447 at 

ProSAPiPl 

ProSAPiPl protein 

0.1 

-0.8 

0 

-1.2 

functions 









204492 at 

ARHGAP11A 

KIAA0013 gene product 

0.5 

1 

-0.2 

0.5 


208154 at 

LOC51336 

Mesenchymal stem cell protein DSCD28 

-0.7 

-0.3 

-0.1 

-1.2 


212686 at 

PPM1H 

Protein phosphatase 1H (PP2C domain 

0 

-0.4 

-0.1 

-1.1 




containing) 






212845 at 

SAMD4 

Sterile alpha motif domain containing 4 

0.1 

-0.4 

1 

-0.5 


213069 at 

HEG'’ 

HEG homolog 

0.3 

-0.3 

1.1 

-0.2 


213256 at 

MGC48332 

Hypothetical protein MGC48332 

0 

0.2 

0.4 

1.1 


214077 x at 

MEIS3 

Meisl, myeloid ecotropic viral integration site 1 

0.5 

0.4 

0.3 

1.2 




homolog 3 (mouse) 






214862 x at 


MRNA; cDNA DKFZp564G1162 

-0.1 

-0.1 

-1 

-0.6 


215281 x at 

POGZ 

Pogo transposable element with ZNF domain 

0 

1 

0.1 

-0.3 


216715 at 



1.2 

1.3 

1 

1.1 


217047 s at 

FAM13A1'’ 

Family with sequence similarity 13, member A1 

-0.3 

0.6 

-1 

-0.1 


217300 at 



-1 

-0.2 

-1 

-0.4 


218031 s at 

C14orfll6 

Chromosome 14 open reading frame 116 

0.1 

-0.7 

0.2 

-1 


218541 s at 

C8orf4 

Chromosome 8 open reading frame 4 

0.4 

0.1 

1.3 

0.9 


218816 at 

LANO 

LAP (leucine-rich repeats and PDZ) and no PDZ 

0.1 

0.3 

0 

1 




protein 






219010 at 

FLJ10901 

Hypothetical protein FLJ10901 

-0.1 

0.1 

1 

1.6 


219166 at 

C14orfl04 

Chromosome 14 open reading frame 104 

-0.4 

0.3 

-1.3 

0.4 


219696 at 

FLJ20054 

Hypothetical protein FLJ20054 

0.1 

-0.1 

-0.2 

-1 


220041 at 

SMP3 

SMP3 mannosyltransferase 

-0.4 

-1.1 

-0.5 

-1 


220441 at 

FLJ13236 

Hypothetical protein FLJ13236 

0 

0.4 

0.5 

1.4 


220770 s at 

LOC63920 

Transposon-derived Buster3 transposase-like 

-0.6 

1.3 

-0.8 

1 


222196 at 

LOC286434 

Hypothetical protein LOC286434 

1.1 

0.8 

0.9 

2.2 


AFFX-r2- 


Human 18S rRNA 6 

-2.1 

1.7 

0.7 

2.6 


Hsl8Sr 








RNA- 








M x at 








AFFX-r2- 


Human 28S rRNA 6 

-1.4 

0.7 

0.5 

1.2 


Hs28Sr 








RNA- 








3_at 








a Minus, decreased gene expression; no minus, increased gene expression; underlining, two-to fourfold-increased or decreased gene expression; boldface, ^-fold- 
increased or decreased gene expression. 

b More than one probe set has a significant change in gene expression; only the most representative is shown. 
c HIV-1, human immunodeficiency virus type 1. 


tried to study the difference in the expression profile at the 
relatively early stages of the infection of 2 and 4 h postinocu¬ 
lation, as reported for herpes simplex virus type 1 (17). This 
time frame is biologically relevant because proliferation of the 
Golgi complex and related vesicles and swelling of trans-Golgi 
sacs were observed in infected cells within the first hour of 
infection. Extracellular virus particles were present in 5% and 
30% of the cell populations at 5 and 6 h postinoculation, 
respectively (35). This would also facilitate the analysis, since a 
lower number of altered genes were involved. In the compar¬ 
ative transcriptomic study, far more genes (n = 136) were 
upregulated by SARS-CoV than by HCoV-229E. In contrast to 
the reported findings of increased antiapoptoptic/inflamma¬ 


tory gene expression and decreased proapoptotic/inflamma- 
tory gene expression in enterocyte cell lines (10), far more 
proapoptotic and proinflammatory genes were expressed in 
Huh7 cells infected by SARS-CoV but not HCoV-229E. For 
instance, expression of BCL2 was induced by SARS-CoV in 
enterocytes, yet we observed upregulation of its antagonists, 
including BAX and BCL2L11, in Huh7 cells. Moreover, much 
higher expression of other proapoptotic proteins, including 
CASP7, CARD 10, PMAIP1, and GADD45B, was also induced 
by SARS-CoV than by HCoV-229E. Furthermore, there was 
marked perturbation of genes involved in cell cycle regulation, 
including induction of the CDKN2B, gene which can mediate 
growth arrest at the G x phase. The induction of proinflamma- 
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FIG. 1. Linkage of the altered genes by using Pathway Assist (version 2.53) shows clustering of genes related to apoptosis, inflammation, and 
coagulation. Only genes with known direct interactions are shown. The software is available from Ariadne Genomics Inc. +, upregulated/increased 
expression; —, downregulated/decreased expression. 


tory cytokines by SARS-CoV was even more prominent when 
compared to HCoV-229E. The induction of IL-8 may be of 
pathogenic importance, as its level has been positively corre¬ 
lated with disease severity in pulmonary infection by respira¬ 
tory syncytial virus. Thus, the observed significantly higher 




TFIP2 

PAI1 

THBS1 

IL8 

NFKB2 

JUNB 

PHLDIA 

CARD10 

BAX 

PBGD 


> 


Coagulation 



Inflammation 

Transcription 

Apoptosis 

House keeping gene 


FIG. 2. Validation of microarray data by semiquantitative RT- 
PCR. The expression patterns of all nine selected genes and the house¬ 
keeping gene, which were assessed for their relative abundances in 
mock infection, at 2 h and 4 h after SARS-CoV inoculation, and at 2 h 
and 4 h after HCoV-229E inoculation, were consistent with the ex¬ 
pression pattern shown in the microarray experiment. 


level of IL-8 induced by SARS-CoV compared with HCoV- 
229E in Huh7 cells may recapitulate the host response to these 
viruses by pneumocytes. The induction of various chemokines 
of the CXC or CCL family may mediate the chemotaxis of 
lymphocytes and neutrophils. These alterations in gene expres¬ 
sion are in keeping with the histological changes of SARS 
hepatitis, in which cellular apoptosis, marked accumulation of 
cells in mitosis with ballooning degeneration of hepatocytes, 
and moderate lymphocytic infiltration were found in liver tis¬ 
sue biopsy samples (4). Recently, a vaccine study in ferrets 
immunized with the modified vaccinia virus Ankara carrying 
the spike protein showed that the ferrets developed hepatitis 
after being challenged with wild SARS-CoV. This finding tends 
to suggest that the Huh7 cell line model may have some rele¬ 
vance to the pathogenesis of SARS (50), but our present find¬ 
ings still may not be directly applicable to the pneumocytes 
involved in SARS. 

The upregulation of genes involved in procoagulation and 
platelet activation was interesting. TFPI2 inhibits thrombin 
generation by binding and inactivation of the tissue factor- 
factor Vila complex. Upregulation of the gene probably rep¬ 
resents an inhibitory response to restrain the activation of the 
coagulation pathway during acute inflammation. In contrast, 
TFPI2 also inhibits both free and matrix/cell-associated plas- 
min, thus favoring fibrin deposition, and may have positive role 
in matrix turnover (11). Upregulation of the PAH gene ac¬ 
companied by a dramatic increase in the protein level results in 
an antifibrinolytic response, favoring fibrin deposition during 
the acute inflammatory phase of the disease. It is important to 
remember that mouse hepatitis virus can activate the immune 
coagulation system by the fgl2 gene, encoding a prothrombi- 
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1000 


♦ SARS-CoV 2 hr ■ SARS-CoV 4 hr A HCoV-229E 2 hr X HCoV-229E 4 hr 



l 


{ 


} 


0.1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 

TFPI2 PAI1 THBS1 It.8 NFKB2 JUNB PHLDIA CARD 10 BAX 


Target gene 


FIG. 3. Fold changes of expression of nine selected genes by RT-qPCR. The P values of the differences in fold changes of the target gene 
expressions are as follows: TFPI2 (SARS-CoV at 2 h versus HCoV-229E at 2 h), 0.02; TFPI2 (SARS-CoV at 4 h versus HCoV-229E at 4 h), 0.02; 
PAI1 (SARS-CoV at 2 h versus HCoV-229E at 2 h), 0.05; PAI1 (SARS-CoV at 4 h versus HCoV-229E at 4 h), 0.05; THBS1 (SARS-CoV at 2 h 
versus HCoV-229E at 2 h), 0.03; THBS1 (SARS-CoV at 4 h versus HCoV-229E at 4 h), 0.35; IL-8 (SARS-CoV at 2 h versus HCoV-229E at 2 h), 
<0.01; IL-8 (SARS-CoV at 4 h versus HCoV-229E at 4 h), <0.01; NFKB2 (SARS-CoV at 2 h versus HCoV-229E at 2 h), 0.05; NFKB2 
(SARS-CoV at 4 h versus HCoV-229E at 4 h), <0.01; JUNB (SARS-CoV at 2 h versus HCoV-229E at 2 h), <0.01; JUNB (SARS-CoV at 4 h versus 
HCoV-229E at 4 h), 0.28; PHLDIA (SARS-CoV at 2 h versus HCoV-229E at 2 h), 0.20; PHLDIA (SARS-CoV at 4 h versus HCoV-229E at 4 h), 
<0.01; CARD 10 (SARS-CoV at 2 h versus HCoV-229E at 2 h), 0.03; CARD10 (SARS-CoV at 4 h versus HCoV-229E at 4 h), 0.04; BAX 
(SARS-CoV at 2 h versus HCoV-229E at 2 h), 0.04; BAX (SARS-CoV at 4 h versus HCoV-229E at 4 h), 0.24. 



—SARS-CoV 2 hr 
—•— SARS-CoV 4 hr 
SARS-CoV 12 hr 
—•— SARS-CoV 24 hr 

- - ♦ - - HCoV-229E 2 hr 

- - • - * HCoV-229E 4 hr 

- - * - • HCo V-229E 12 hr 

- - ♦ - - HCoV-229E 24 hr 


FIG. 4. Validation of protein expression for TFP1, PAI1, and IL-8 by enzyme immunoassay. The P values of the differences in protein expression are 
as follows: TFPI2 (SARS-CoV at 2 h versus HCoV-229E at 2 h), <0.01; TFPI2 (SARS-CoV at 4 h versus HCoV-229E at 4 h), 0.37; TFPI2 (SARS-CoV 
at 12 h versus HCoV-229E at 12 h), 0.18; TFPI2 (SARS-CoV at 24 h versus HCoV-229E at 24 h), 0.12; PAI1 (SARS-CoV at 2 h versus HCoV-229E 
at 2 h), <0.01; PAI1 (SARS-CoV at 4 h versus HCoV-229E at 4 h), <0.01; PAH (SARS-CoV at 12 h versus HCoV-229E at 12 h), <0.01; PAH 
(SARS-CoV 24 at h versus HCoV-229E at 24 h), 0.01; IL-8 (SARS-CoV at 2 h versus HCoV-229E at 2 h), <0.01; IL-8 (SARS-CoV at 4 h versus 
HCoV-229E at 4 h), <0.01; IL-8 (SARS-CoV at 12 h versus HCoV-229E at 12 h), <0.01; IL-8 (SARS-CoV at 24 h versus HCoV-229E at 24 h), <0.01). 
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nase (12). This enzyme can induce macrophage procoagulation 
activity which results in fibrin deposition on the endothelia of 
intrahepatic veins and hepatic sinusoids. The end result could 
be confluent hepatocellular necrosis. The low number of liver 
biopsies done with these patients may account for the lack of 
reports on these full-blown changes related to vascular dam¬ 
age. However, systemic vasculitis, including edema, localized 
fibrinoid necrosis, and infiltration by monocytes, lymphocytes, 
and plasma cells into vessel walls in various tissues, was re¬ 
ported (13). Thrombosis was found in small veins. Interest¬ 
ingly, marked upregulation of a proapoptotic gene, PHLDA1, 
was observed in SARS-CoV infection of Huh7. A previous 
study has shown that overexpression of this gene in vascular 
endothelial cells would lead to decreased cell adhesion and 
induce detachment-mediated apoptosis. This gene, if similarly 
induced in vascular endothelial cells infected by SARS-CoV, 
may contribute to the vascular damage induced by SARS-CoV 
infection (21). 

The clinical manifestations and histological changes involv¬ 
ing the pneumocytes, enterocytes, and vascular endothelium 
were not surprising because ACE2, the host cell receptor for 
viral entry, was present in all these cell types (19). Similarly, 
abnormal urinalysis associated with increased viral load in 
urine is expected, since SARS-CoV could be cultured from or 
antigen expression could be detected in kidney cell lines and 
kidney tissues of patients (22). However, neither ACE2 nor 
L-SIGN could be detected by immunohistochemical staining in 
normal human hepatocytes or Huh7 (1, 42). This could be 
explained by an alternative cellular receptor for viral entry or 
just because the amount of receptor required for cell entry is so 
low that the cells may appear negative on immunohistochem¬ 
ical staining. Our previous study has shown that a high viral 
load in serum or stool is associated with hepatic dysfunction 
(22). The viral load in serum reflects the dynamic of the viral 
replication from any organs and the clearance mechanism of 
the host, whereas a high viral load in the stool may indirectly 
reflect the degree of portal venous viremia. Thus, the lack of 
cell death and inflammation in the intestine may still have 
importance in the pathogenesis of SARS in other body sites. 

In conclusion, SARS-CoV produces more severe distur¬ 
bance of host cell gene expression in a human epithelial cell 
line of liver origin than does HCoV-229E at the early stage of 
infection. There are marked alterations in gene expression 
related to apoptosis, inflammation, and procoagulation. These 
findings are consistent with the histological changes of SARS, 
especially in the liver and blood vessels. Besides antivirals (6, 
24), other modalities of treatment, such as antiapoptotic agents 
(34), immunomodulators against inflammation (7), and modi¬ 
fiers of coagulation, should be considered in future research. It 
is important to remember that many patients continued to 
deteriorate despite a decreasing viral load 2 to 3 weeks after 
the onset of SARS. 
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